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ABSTRACT 

Context. Tight binaries discovered in young, nearby associations, with known distances, are ideal targets to provide dynamical mass 
measurements to test the physics of evolutionary models at young ages and very low masses. 

Aims. We report for the first time the binarity of TWA22. We aim at monitoring the orbit of this young and tight system to determine 
its total dynamical mass using an accurate distance determination. We also intend to characterize the physical properties (luminosity, 
effective temperature and surface gravity) of each component based on near-infrared photometric and spectroscopic observations. 
Methods. We use the adaptive optics assisted imager NACO to resolve the components, to monitor the complete orbit and to obtain 
the relative near infrared photometry of TWA22 AB. The adaptive optics assisted integral field spectrometer SINFONI was also used 
to obtain medium resolution {R^ = 1500 - 2000) spectra in JHK bands. Comparison with empirical and synthetic librairies were 
necessary to derive the spectral type, the effective temperature and the surface gravity for each component of the system. 
Results. Based on an accurate trigonometric distance (17.53 ± 0.21 pc) determination, we infer a total dynamical mass of 220 ± 
21 Mj,ip for the system. From the complete set of spectra, we find an effective temperature T^ff = 2900^,[jQ K for TWA22 A and 
Tfff = 2900^^1515 ^ f"'- TWA22 B and surface gravities between 4.0 and 5.5 dex. From our photometry and a M6 ± 1 spectral 
type for both components, we find luminosities of log(L/Lg)=-2.1 1± 0.13 dex and log(L/LQ)=-2.30± 0.16 dex for TWA22 A and B 
respectively. By comparing these parameters with evolutionary models, we question the age and the multiplicity of this system. We 
also discuss a possible underestimation of the mass predicted by evolutionary models for young stars close to the substellar boundary. 

Key words. Stars: fundamental parameters, low-mass, brown dwarfs, binary (TWA22 AB): close, formation - Instrumentation: adap- 
tive optics, spectrographs 



1. Introduction 

Mass and age are fundamental parameters of stars and brown 
dwarfs that determine their luminosity, effective tempera- 
ture, atmospheric composition and surface gravity commonly 
derived through photometric and spectroscopic observations. 
Evolutionary models are cuiTently widely used in the commu- 
nity to infer masses of stars and brown dwarfs, but they rely 
on equations of states and atmospheric models non calibrated at 
young ages and at very low masses. However, direct mass mea- 
surements can be obtained by the mean of different observing 
techniques, combining light curve studies to radial velocity on 
ecUpsing binaries, astrometric follow-up with double lined spec- 
troscopy of tight binaries or measuring the Keplerian motion of 
circumstellar disks. In recent years, direct mass measurements 
for 23 pre-main sequence stars with masses ranging from 0.5 to 2 



* Based on service-mode observations (072.C-0644, 073.C-0469, 
075.C-0521, 076.C-0554, 078.C-0510, 080.C-058I) collected at the 
European Organisation for Astronomical Research in the Southern 
Hemisphere, Chile. 



Mq showed discrepancies with predictions by up to a factor of 2 
in mass and 10 in ages (Mathieu et al. 2007 1. Such measurements 
are more rare for lower masses (< O.SMq) systems. [Hillenbrand] 
& White (2004 ) showed that the models tend to understimate the 
mass of the companion UZ Tau Eb (M= 0.294±0.027 Mp, age 
'5 Myr, |Prato et al.| ( |2002| )). [Close et aL] ([2005 ) deiived similar 
conclusions but the age and the luminosity of the companion AB 
Do r C (M^0. 090±0.005 Mq, age-75 Myr) is still under debate 
(Boccaletti et al. 2008 ). And recently, the surprising discovery 



of the unpredicted temperature reversal ( Stassun et al.|2007] l be 
tween 2M035 A (M=0.0541 ±0.0046 Mq, age --1 Myr) and its 
companion (M=0.0340±0.0027 Mq, age -1 Myr) proves the ne- 
cessity to find more calibrators. The challenge is to determine 
unambiguously their physical properties (mass, L, T^jj, g and 
age) and to explore as much as possible the parameter space cov- 
ered by evolutionary models. The influence of other parameters 
such as metallicity needs also investigation ( Boden et al.|2005" ^ 
Burgassei',2007j . 
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Table 1. Observing log. 



UT Date 


Name 


Instrument 


Mode 


Filter 
(Grism) 


Camera 


Airmass 


Seeing 
(arcsec) 


EC|^ 

(%) 


Exp. 
Time 

(s) 


Note 


2004/03/05 


TWA22AB 


NACO 


imaging 


NB^.17 


S27 


1.16 


1.23 


16.5 


20 




2004/03/05 


TWA22AB 


NACO 


imaging 


NBJ.24 


S13 


1.15 


1.40 


27.7 


120 




2004/04/27 


TWA22AB 


NACO 


imaging 


NBJ.75 


827 


1.16 


0.84 


48.2 


6 




2005/05/06 


TWA22AB 


NACO 


imaging 


H 


S13 


1.15 


0.72 


12.5 


5 




2006/01/08 


TWA22AB 


NACO 


imaging 


J 


S13 


1.16 


0.84 


17.5 


50 




2006/01/08 


TWA22AB 


NACO 


imaging 


H 


S13 


1.17 


0.65 


17.4 


50 




2006/02/26 


TWA22AB 


NACO 


imaging 


H 


S13 


1.15 


1.10 


32.1 


150 




2006/02/26 


TWA22AB 


NACO 


imaging 


Ks 


S27 


1.15 


1.50 


10.2 


50 




2007/03/06 


TWA22AB 


NACO 


imaging 


H 


S13 


1.15 


0.60 


45.7 


103.5 




2007/12/04 


TWA22AB 


NACO 


imaging 


H 


S13 


1.28 


1.00 


34.4 


150 




2007/12/26 


TWA22AB 


NACO 


imaging 


H 


S13 


1.16 


0.90 


39.9 


150 




2007/02/12 


TWA22AB 


SINFONI 


spectroscopy 


J (2000) 


25 


1.18 


1.56 


12.1 


1080 




2007/02/13 


TWA22AB 


SINFONI 


spectroscopy 


J (2000) 


25 


1.16 


0.98 


18.5 


1080 




2007/02/13 


HIP049201 


SINFONI 


spectroscopy 


J (2000) 


25 


1.13 


0.86 


27.0 


20 


Telluric Standard 


2007/02/13 


HIP038858 


SINFONI 


spectroscopy 


J (2000) 


25 


1.14 


0.78 


28.5 


20 


Telluric Standard 


2007/02/12 


HIP035208 


SINFONI 


spectroscopy 


J (2000) 


25 


1.17 


1.17 


15.2 


60 


Telluric Standard 


2007/02/11 


TWA22AB 


SINFONI 


spectroscopy 


H+K (1500) 


25 


1.16 


0.89 


29.5 


960 




2007/02/09 


TWA22AB 


SINFONI 


spectroscopy 


H+K (1500) 


25 


1.52 


1.00 


17.2 


960 




2007/02/11 


HIP052202 


SINFONI 


spectroscopy 


H+K (1500) 


25 


1.19 


0.93 


24.7 


20 


Telluric Standard 


2007/02/09 


HIP052202 


SINFONI 


spectroscopy 


H+K (1500) 


25 


1.55 


1.22 


14.0 


40 


Telluric Standard 



" Corresponds to the mean strehl ratio for the spectroscopic observations. 



The TW Hydrae association (TWA) is the first co-moving 
group of young (< 100 Myr), nearby (< 100 pc) stars, that 
was identified near the Sun ( |Kastner et al.|[T997] l. Ideal obser- 
vational niche for the study of stellar and planetary formation, 
TWA was actually the tip of an iceberg composed of hundreds of 
young stars, spread in different groups, that were identified dur- 
ing the last decade ( [Zuckerman & Soiig 2004 Torres et al. 2008| l. 
TWA counts now 27 members covering a mass regime from 



intermediate-mass stars to planetary mass objects (Chauvin et al. 
2005| l. Its 8.3 + 0.8 Myr dynamical age was found by a conver- 
gence method (|de la Reza et al.|[2006| . Independently, Barrado 



phot ' 

poo 



|Y Navascues (2006 1 estimated an age of 10^^ Myr from the 
phot ometry, the activity and the lithium depletion. IScholz et al. 



2007 ) show that the association is Myr old by comparing 
rotational velocities with published rotation periods for a subset 
of stars. Finally, |Mentuch et al. ( 20081 estimated an age of 12 ± 
8 Myr from the study of lithium depletion in five nearby young 
associations (hereafter M08). 

Song et al.| ( |2003) , hereafter S03, identified TWA22 as an 



M5 member of TWA. The strong Li /16708 A feature supported 
the extreme youth of this member Later, Mamajek'(2005) ques- 
tioned the membership of TWA22 from a kinematic study of 
TW Hydrae members. Finally, [Song et aL) ( |2006| l discussed the 
Mamajek (2005j) results that appeared in disagreement with the 
very strong lithium line of the source. The proximity (17.5 + 0.2 
pc, seeTeixe ira et al.|p009| ) and the reported youth of TWA22 
by S03 made it consequently a potential target for the detection 
of companions at small orbital radii. 

In our program for detection of companions in young associ- 
ations using the adaptive optic (AO) assisted imager NACO, we 
resolved TWA22 as a tight (~100 mas) binary. With a projected 
physical separation of 1.76 ± 0.10 AU (see Fig.[T]i, this system 
offered a unique opportunity to measure its dynamical mass and 
to possibly test the evolutionary model predictions at young ages 
using combined photometric and spectroscopic observations. 

We report here the discovery of the TWA22 binarity and the 
results of a dedicated 4 years observing program, using com- 




Fig. 1. VLT/NACO image of TWA22 AB obtained in H-band 
with the S13 camera on December 26th, 2007. North is up and 
East is left. The field of view is 1 " x 1 ". 



bined imaging and 3D-spectroscopy with AO. The purpose was 
to measure the dynamical mass of TWA22 AB and to charac- 
terize the physical properties of the individual components. In 
§2, we describe our AO observations with the VLT/NACO im- 
ager and with the VLT/SINFONI integral field spectrograph. The 
associated data reduction and spectral extraction techniques are 
detailed in §3. In §4, we present our orbital solutions and our 
spectral analysis. In §5 we compare and discuss the evolution- 
ary model predictions associated to our dynamical mass mea- 
surement with the physical properties (surface gravity, temper- 
ature and luminosity) derived from our photometric and spec- 
troscopic observations. This leads us to discuss the membership 
of TWA22 AB to the TW Hydrae association, the multiplicty of 
the system and a possible underestimation of the mass predicted 
by evolutionary models for young stars close to the substellar 
boundary. 
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2. Observations 

2.1. VLT/NACO Observations 

TWA22 AB was observed at the 8.2 m VLT UT4 Ye pun with 
the Nasmyt h Adaptive Optic (AO) system NAOS ( |Rousset 
|et al.|| 2000| coupled to the High-Resolution Near-IR Camera 
CONICX pIenzen et al.|[T998l ). NAOS and CONICA (NACO) 
resolved the system as a tight binary for the first time on March 
5th, 2004. Follow-up observations were conducted during four 
years from early 2004 to end 2007. To image TWA22AB, 
we used the narrow band filters: NBJ.24 {Ao-l.231 fim, 
FWHM=0.015 yum), NBJ.75 (^o=l-748 jim, FWHM=0.026 
A<m), NB^.17 (/lo=2.166 ^im, FWHM=0.023 yum). The broad 
band filters J {Aq=\.21 fim, FWHM=0.25 fim), H (^o=l-66 jjm, 
FWHM=0.33 jum) and K, (/lo=2.18 fim, FWHM=0.35 fim) were 
also used coupled to a neutral density (attenuation factor of 80). 
CONICA was used with the S13 and S27 cameras to Nyquist- 
sample the PSF depending on the selected filter. The data were 
recorded under seeing ranging from 0.6" to 1.5" (see Table [TJ. 
TWA22 AB was bright enough in visible to be used by NAOS 
for wave-front analysis. For each observation period, dithering 
around the object in J, H and Ks bands combined with nodding 
were necessary to run a good sky estimation during the data re- 
duction process (see part 2.2). PSF references were observed at 
different airmasses with identical setups. The ©i Ori C astromet- 
ric field (McCaughrean & Stauffer 1994 ) was observed at each 
epoch to calibrate the detector platescale and orientation when- 
ever necessary. The results are reported in Table |2] 



Table 2. Mean plate scale and detector orientation for our differ- 
ent observing NACO runs. 



UT Date 


Camera 


Plate Scale 


Orientation of true north 








east of the vertical 






(mas/pixel) 


(0) 


2004/03/05 


S27 


27.01 ± 0.05 


-0.18 ±0.19 


2004/04/27 


S27 


27.01 ±0.10 


0.08 ± 0.20 


2005/05/06 


S13 


13.25 ±0.08 


-0.02 ±0.10 


2006/01/08 


S13 


13.26 ±0.08 


0.18 ±0.10 


2006/02/26 


S13 


13.25 ± 0.06 


0.19±0.11 


2007/03/06 


S13 


13.28 ±0.10 


0.08 ±0.15 


2007/12/04 


S13 


13.29 ±0.12 


-0.21 ±0.17 


2007/12/26 


S13 


13.29 ±0.12 


-0.21 ±0.17 



2.2. VLT/SINFONI observations 

The SINFONI instrument (Spectrograph for INtegral Field 
Observations in the Near Infrared), located at the Cassegrain fo- 
cus of the VLT UT4 Yepun, was used to observe TWA22AB 
between February 9 and 13 2007. SINFONI includes an integral 
field spectrometer SPIFFI (Spectrograph for Infrared Faint Field 
Imaging, see jEisenhauer et al.| ^2003)), operating in the near- 
infrared (1.1 - 2.45 fim). SPIFFI is assisted with the 60 actuators 
Mutlti-Applications Curvature Adaptive Optic system MACAO 
([Bonnet et alT 2003). We used the small SPIFFI field of view of 
0.8" X 0.8" corresponding to a plate scale of 25 mas per pixel 
to Nyquist-sample the SINFONI AO coiTected PSF. The field of 
view is optically sliced into 32 horizontal slitlets, that sample the 
horizontal spatial direction and which are reaiTanged to form a 
pseudo-long slit. Once dispersed by the grating on the 2048 x 
2048 SPIFFI detector, each shtlet of 64 pixels width (spatial di- 
rection) corresponds to 64 spectra of 2048 pixels long (spectral 



direction). The 2048 independent spectra on the detector are re- 
organized during the reduction process in a datacube which con- 
tains the spatial (X, Y) and the spectral (Z) informations. The 
cube is resampled in the vertical dimension (Y) to have the same 
number of pixels as in X. 

To cover the full spectral range between 1.1 to 2.45 fim, in- 
dividual integrations times of 90 s were necessary to image the 
system in the J band (1.1 - 1.4 fim, R=2000) and 20 s in the 
H+K band (1.45 - 2.45 fim, R=1500). For each band, dithering 
around the object was used to increase the field of view and to 
supress residual bad pixels, leading to a total observing time on 
target of ~ 5 min. An additional frame was acquired on the sky to 
improve our correction. The adaptive optic loop was locked on 
TWA22 AB itself Standard stars HIP038858 (B3V), HIP049201 
(B2V), HIP035208 (B3 V) and HIP052202 (B4V) were observed 
at similar airmasses to remove the telluric lines (see Table [TJ- 



3. Data reduction and analysis 

3. 1 . Higii Contrast Imaging 

For each observation periods, the ESO eclipse reduction soft- 
ware (Devillard 1997 1 dedicated to AO image processing was 
used on the complete set of raw images. Eclipse computes bad- 
pixel detection and interpolation, flat field correction and aver- 
aging pairs of shifted images with sub-pixel accuracy. The soft- 
ware run sky estimation on object-dithered frames using median 
filtering through the frame sequence. 

A deconvolution algorithm dedicated to stellar field blurred 
by the adaptive optics coiTected point spread functions ( Veran & 
|Rigaut'1998) was applied on TWA22 AB images to accurately 
find the position and the photometry of the companion relative 
to the primary. The algorithm is based on the minimization in 
the Fourier domain of a regularized least square objective func- 
tion using the Levenberg-Marquardt method. We used Nyquist- 
sampled unsaturated images of standard stars obtained the same 
night as TWA22 observations with identical setups under vari- 
ous atmospheric conditions. These frames captured the variation 
of AO corrections. They were used as input point spread func- 
tions (PSF) to estimate the decon volution process erro r The IDL 
Starfinde^PSF fitting package ( Diolaiti et al.|2000 i confirmed 
these results. 



IDL procedures can be downloaded 
|http://www.bo. astro.it/~giangi/StarFinder/index.htm 



at 



Table 3. Relative positions and contrasts of TWA 22 A and B. 
Magnitude difi'erences are given in the NACO photometric sys- 
tem. 



UT Date 


Filter 


Camera 


^0-72000 


A(5y2000 


Am 










(mas) 


(mas) 


(mag) 




2004/03/05 


NB2.17 


S27 


99 ±3 


-17 ±3 


0.33 ± 


0.12 




NB1.24 


S13 






0.37 ± 


0.10 


2004/04/27 


NB1.75 


S13 


98 ±6 


-36 ±6 


0.41 ± 


0.10 


2005/05/06 


H-ND 


S13 


15 ±3 


-89 ±3 


0.66 ± 


0.05 


2006/01/08 


J 


S13 






0.40 ± 


0.04 




H 


S13 


-68 ±2 


-49 ±2 


0.54 ± 


0.05 


2006/02/26 


H 


S13 


-74 ± 3 


-30 ±3 


0.54 ± 


0.05 




K, 


S27 






0.46 ± 


0.18 


2007/03/06 


H 


S13 


-57 ±4 


80 ±2 


0.49 ± 


0.03 


2007/12/04 


H 


S13 


19 ±3 


98 ±3 


0.52 ± 


0.10 


2007/12/26 


H 


S13 


26 ±3 


97 ±3 


0.53 ± 


0.10 
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3.2. Integral Field Spectroscopy 

We used the SINFONI data reduction pipeline (1 .7. 1 version, see 
|Modigliani et al.| ( |2007 !)) for raw data processing. The pipeline 
carries out cube reconstruction from raw detector images. The 
flagging of hot and non-linear pixels is executed in a similar way 
as in NACO images. The distortion and wavelength scale are 
calibrated on the entire detector using arc-lamp frames. Slitlet 
distances are accurately measured with North-South scanning 
of the detector illuminated with an optical fiber. In the case 
of standard stars observation, object-sky frame pairs are sub- 
tracted, flat fielded and corrected from bad pixels and distor- 
tions. Datacubes are finally reconstructed from clean science im- 
ages and are merged in a master cube. Spectra of standard stars 
cleaned from stellar lines are finally used to correct the TWA22 
AB spectra from telluric absorptions. 

TWA22 A and B are centered and oriented horizontally in 
the J and H+K master cubes with a field of view of 1 . 1 " x 1.1". 
Atmospheric refraction induces diff'erent sources positions for 
different wavelengths within the instrument field of view and in- 
creases with airmass. Combined with the small SINFONI field 
of view, this produces differential flux losses that were noticed 
in the bright standard stars datacubes. This effect remains lim- 
ited for TWA22. The cubes of 1 1 February 2007 appear to have 
some spaxels contaminated by flux oscillations of a few ADUs. 
These oscillations were not negligible and blurred CO bands at 
2.3 yum. They are present along the dispersion axis in the raw de- 
tector images of both HIP052202 and TWA22. Their amplitudes 
do not remain constant in time but follow a 15.3 pixels period. 
We then filtered partially this contribution on each individual im- 
age in the Fourier space using a pass-band function. The origin 
of the problem is likely to be related to 50 Hz pick-up noise. 

We used a modified version (|Dumas et al.||2001) of the 
CLEAN algorithm (Hogbom 1974; Schwartz 1978) to extract 
separately the flux of TWA22 A and B in each monochromatic 
images contained in the datacubes. The standard star are used 
as initial PSF-references. Once scaled to match the TWA22 A 
maximum at the primary position and for all wavelengths, the 
PSF is subtracted to the TWA22 AB datacube. The sequence is 
repeated to model the secondary contribution, cleaned from the 
primary wings, and to provide a new PSF-reference. After a few 
iterations minimizing the final quadratic residual datacube, the 
spectra of each individual component are extracted. 

The algorithm was first adapted to work on cube images. 
Unfortunately, the difference of sampling between the X and 
Y directions limited the sub-pixel shift accuracy. We therefore 
collapsed the cube along the Y-axis in order to obtain the flux 
profile along the X direction. We chose to duplicate the primary 
flux profile for the PSF model. The algorithm converged in a 
few iterations and produced extracted spectra in J and H+K with 
an extraction error less than 5%. The extracted spectra were di- 
vided by standard star spectra corrected from intrinsic features 
and multiplied by a black body spectrum at the standard star 
temperature. The SINFONI pipeline coefficients were used for 
wavelength calibration. 

4. Results 

4.1. Astrometry, Orbit and dynamical mass 

The relative positions of TWA 22 A and B (B with respect to 
A) at all observation epochs are reported in Table [3] The data 
allows a determination of the mutual orbit of the binary. We de- 
fine a cartesian referential frame (O, X, Y, Z) where X points to- 
wards the north, Y toward the east and Z toward the Earth. The 



(OXY) plane corresponds thus to the plane of the sky. Then in a 
Keplerian formalism, the (x = A6,y = Aa) projected position of 
the binary onto the plane of the sky reads 

a 

X - — (cos u - e) (cos oj cos Q - sin w cos / sin Q) 
d 



H — Vl - sin M (- sin w cos Q. - cos o) cos / sin Q) , (1) 
d 
a 

y = - (cos u - e) (cos w sin Q + sin w cos / cos Q.) 
d 



H — Vl - e^smu (- sin sin O + cos oj cos Q. cos /) , (2) 
d 

where a is the semi-major axis of the orbit (in AU), d is the 
distance of the binary (in pc), e is the eccentricity, / is the incli- 
nation, Q is the longitude of ascending node (counted from north 
towards east), u) is the argument of periastron, and u is the ec- 
centric anomaly that describes the current location of the binary 
along its orbit, u is related to the time t by the classical Kepler's 
equation 



2;r 



-{t - t^) — u — e sin u 



(3) 



where T is the orbital period and fp is the time reference for peri- 
astron passage. Once the distance d is known (17.5 pc, Teixeira 
|et al.| ( |2009[ l), the fit of the observational data allows to deter- 
mine the 7 parameters T, a, e, i, Q, co and fp. Then Kepler's third 
law leads to the determination of the total mass M. 

The fit is performed via a Levenberg-Marquardt;^^^ minimiz- 
ing algorithm. In practice, instead of (e, /, Q., a>), the equations 
are solved for the classical variables 

/ 

k — e cos(f2 + oj) q - sin - cos Q 

/ 

h - e sin(f2 + tj) p - sin - sin f2 . (4) 

which avoids singularities towards small eccentricities and incli- 
nations. The uncertainties on the fitted parameters are estimated 
from the resulting covariance matrix at the end of the fit proce- 
dure. 

Levenberg-Marquardt is an interative gradient method for 
converging towards a mininum of the function. Depending 
on the starting guess point, many local minima can be found. In 
the present case, all the attemps we made (by letting the start- 
ing point vary) converge towards the same solution that is listed 
in Table |4] and viewed in projection onto the plane of the sky in 
Fig.|2] The available astrometric data set appears to cover almost 
one complete orbital period with a good sampling of the perias- 
tron passage. We are thus confident in our fitted solution. The 
orbit appears then slighly eccentric (e ^ 0.1) and viewed close 
to pole-on / ^ 27° from the Earth. 



4.2. Photometry 

Table |3] summarizes the magnitude differences between TWA22 
A and B, measured with NACO at different wavelengths. 
Taking into account the filter tranformations between NACO and 
2MASS and the photometry of the unresolved system given from 
the 2MASS Survey ( jCutri et aL|2 003), we derived the apparent 
JHK magnitudes of each component (see Table|5]l. Observations 
under bad seeing conditions were excluded. Based on an accu- 
rate distance (17.53 ± 0.21 pc) determination (Teixeir a et al.| 
2009 ), the absolute magnitudes were also derived. TWA22 be- 
ing a young mid-M system, we monitored its photometric varia- 
tions in the H band. We only noticed a 0.05 variations of the total 
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100 



-100 



a — 


1.77 AU 




d = 


-- 17.53 pc 


e = 


0.099 




M = 


219.49 Mj 


i = 


27.43° 




T = 


= 5.144 yr 


- Q - 








,^ — ^ ono 




134.95° 26/i2/2oo;yii*'*7T 






CJ = 


99.31° 




1^06/03/2007 


tp = 


2006. 03a^ 










Jff05/03/E004 










k_27/04/2004 




5f 26/02/2006 








ife08/01/2006 




1 . . . . 


1755/2005 

1 







Table 5. TWA22 A and B individual magnitudes converted into 
the 2MASS system. 



too 



Aa (mas) 



-100 



Fig. 2. Orbital fit of the relative positions of TWA22 AB ob- 
served from March 2004 to December 2007, as projected onto 
the plane of the sky. The crosses represent the observational data 
with their error bars, the solid line is the fitted projected orbit, 
and the dots correspond to the predicted positions of the model 
at the times of the observations. The dashed line sketches the 
projected direction of the periastron of the orbit. On Jan. 8, 2006, 
the binary was actually very close to periastron. 

Table 4. TWA22 B orbital parameters as determined from the fit 
of the astrometric data (see text for the definition of the parame- 
ters). The reduced of the fit is also reported. 



Reduced 


1.127 


a(AU) 


1.77 ±0.04 


e 


9.95x10-2 ±0.32x10-2 


i (degrees) 


27.43 ± 4.40 


CI (degrees) 


134.95 ± 0.40 


(i> (degrees) 


100 ± 10 


Tp (years) 


2006.039 ± 0.010 


T (years) 


5.15 ±0.09 


M {Mj„„) 


220 ± 21 



magnitude of the system along time. This variation is reported in 
the error bars on our photometry in Table |5] 

Reported in a color (J - K) -magnitude (Mxs) diagram, the 
TWA22A and B photometry can be compared with the pho- 
tometry of M dwarfs of the young, nearby associations TW 
Hydrae (~ 8 Myr), /3 Pictoris (~ 12 Myr), Tucana-Horologium 
(~ 30 Myr) and AB Doradus (~ 70 Myr). Predictions of evolu- 
tionary models of ( |Baraffe et al.|1998| also named NEXTGEN) 
are also given at these young ages (see Fig. |3]l. Although age- 
dependent, the near-infrared photometry of TWA22 A and B ap- 
pears compatible for both components with that expected for 
young mid-M dwarfs but does not allow to give an age estima- 
tion for the binary. The NEXTGEN tracks also appears bluer 
than 10 Myrs old mid-M dwarfs by ~ 0.2 mag, which might be 
related to a non full representation of their spectral energy distri- 
bution. We provide in the following an improved estimation of 
the spectral type of our targets, using our spectroscopic data. 
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8.61 ± 


0.15 
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0.15 


7.39 ± 


0.18 


7.90 ± 


0.18 


Ks 


8.20 ± 


0.19 


8.70 ± 


0.25 


7.02 ± 


0.23 


7.48 ± 


0.28 



TW Hydrae M dwarfs 
p Pic M dwarfs 
Tuc Hor U dwarfs 
AB Doradus M dwarfs 
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0.8 1.0 
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Fig. 3. Color (J - K) - magnitude (Mks) diagram of TWA22 A 
and B compared with the photometry of young M dwarfs mem- 
bers of the TW Hydrae (age=8 Myr, filled triangles), /3 Pictoris 
(age=10 Myr, diamonds), Tuc -Hor (age=30 Myr, stars) and AB 
Dor (age=70 Myr, crosses) associations. Distances and spec- 
tral types were taken from [Torres et al.|p008|l exepted for the 
members of the TW Hydrae moving group dMamajek 2005 i. 
Typical error bars ar e 0.03 for Mk and 0.08 for J-K,. NEXTGEN 
( |Baraffe et aO TW) 8 Myr (dashed-dotted black Une), 10 Myr 
(long-dashed-dotted black line), 30 Myr (long-dashed black line) 
and 70 Myr isochrones (full black line) are overplotted. DUSTY 
tracks ( jChabrier et al.|[2000 ) at identical ages are added (red 
lines). These tracks were shifted by 0.2 mag to redder J-Ks. This 
artificially compensates the remaining incompleteness, at rela- 
tively warm temperatures, in the AMES linelists used in DUSTY 
atmospheric models ( [Mohanty et~aL][2007) l. The NEXTGEN 
tracks also appears bluer than 10 Myrs old mid-M dwarfs by 
~ 0.2 mag, which might be related to a non full representation 
of their spectral energy distribution. At a given age, early-type 
M dwarfs (violet to green points) are bluer and more luminous 
than late-type M dwarfs (green to red points) in the K band. At a 
given spectral type (i.e. temperature) the objects go fainter with 
age. The photometry is compatible with the one of young mid-M 
dwarfs. The unresolved 2MASS photometry is plotted for com- 
parison (black dot). 



4.3. Spectroscopic Analysis 
4.3.1. Line identification 

To identify the numerous spectral features in the TWA22 A and 
B spectra between 1.10 to 2.45 yum, the spectra were compared 
with an homogeneous medium resolution (R^ ~ 2000) sequence 
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of field dwarfs from Gushing et al.| ( [2005] l (hereafter COS; see 
Fig.|4j|5]and|6|l. The TWA22 A and B spectra appear very simi- 
lar. 

In y-band, both TWA22 A and B spectra are dominated by 
the strong Na I doublet at 1.138 fim, the deep K I lines at 1.169, 
1.177, 1.243 and 1.253 fim and the presence of a broad H2O ab- 
sorption from 1.32 to 1.35 fim. Fe I absorptions are also detected 
around 1.170 fim. One is blended with the 1.177 yum K I line. 
We notice additional broad FeH absorptions around 1 .20 yum and 
1 .24 /xm compatible with that expected for mid-M dwarfs as well 
as the presence of the weak Q-branch at 1 .22 fim. Finally, the Al 
I doublet at 1.313 fim is detected. This doublet is expected to 
disappear at the M-L transition. 

In //-band, the spectra are affected by H2O absorptions from 
1.45 to 1.52 fim and from 1.75 to 1.8 fim. They exhibit pro- 
nounced K I atomic lines at 1.517 fxm as well as weak doublets 
of Mg I at 1.503 fim and Al I at 1.675 fim. Weak FeH absorptions 
are also present. They increase from M5 to the M-L transition 
Gushing et al. ( 2003| l and their depths are here compatible with 
those expected for M5 to M7 field dwarfs. 

In /T-band, strong H2O absorptions appear from 1.95 to 2.04 
fim and from 2.3 to 2.45 fim. They are typical from mid-M to 
mid-L dwarfs. Strong Ga I features are present from 1.9 to 2.0 
fim. They tend to disappear in the spectra of field dwarfs at the 
M-L transition. We identify firmly the first overtone of GO near 
2.3 fj.m, the rest being affected by the 50 Hz pick-up noise oscil- 
lations mentionned earlier Additional weak Mn I, Ti I, Mg I and 
Si I absorptions are spread over the J, H and K bands. These lines 
are expected to be rapidly replaced by molecular absorptions for 
dwarfs later than M5. The 1.106 fim band seems to be overlap- 
ping H2O and TiO absorptions with increasing depths from early 
to late M dwarfs. Finally, the 1.626 fim feature corresponds to 
close OH lines, as noted in Leggett et al.| ( |1996] l. 

To conclude, the features detected over the spectra of 
TWA22A and B between 1.1 and 2.45 fim suggest that both 
components have a cool atmosphere, typical of mid to late-M 
dwarfs. 
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Fig. 4. TWA22 A and B spectra compared to spectra of M5 
to M8 dwarfs. The M6 dwarf spectrum reproduces well the J 
bands of TWA22 A and B. However, our spectra seem to have 
a slightly redder slope. We reported identified atomic features 
(blue). Molecular absorptions (FeH bands were identified by 



Gushing et al. (2003 1) are flagged in green and telluric residu- 



als in orange. Atomic absorptions are indicated in blue. 



4.3.2. Continuum fitting and spectral indexes 

The continuum of both TWA22 A and B spectra were compared 
to spectra of field M dwarfs obtained by G05 and [McLean et al. 
( |2003| l, hereafter ML03. Least squares were computed on parts 
of the spectra free from telluric correction residuals. From 1.10 
fim to 1 .27 fim, the TWA22 A and B continuums are well repro- 
duced by M6 + 1 dwarfs. The //-band spectra are visually poorly 
reproduced. Least squares are minimized for M9 dwarfs but with 
2 subclasses uncertainty. Finally, our /T-band spectra are well fit- 
ted by M5 to M7 dwarfs. From these comparisons, we assign a 
spectral type M6 ± 1 to both TWA22 AB components. 

TWA22AB being a young system, we have tested if us- 
ing high surface gravity spectra of old field dwarfs might affect 
our spectral analysis. Intermediate surface gravity reduces the 
strength of alkali lines ( Lucas et al.||2001[ Gorlova et al.|20"0T 



[McGovern et al.|2004[|Kirkpatrick et al.|2606| ) and produces tri- 
angular shape in //-band interpreted as collision induced absorp- 
tions (GIA) of H2. Our spectra were then compared with young 
(age < 8 Myr) dwarfs spectra ( [Gorlova etal . 2003 ; Slesnick et al. 
2^[Lodieu et al.[2008| l in the J and K bands (see Fig.|7]and[8]l. 
They are mostly similar to M5, M6 and M7 dwarf spectra, and 
consistent with the continuum fit obtained with field dwarfs. In 
both cases, our J-band spectra are slightly redder and our H-band 
spectra visually are still poorly reproduced by young and old M 
dwarfs. 



To complete this spectral type determination, spectral in- 
dexes developed by ML03 (from H2O bands at 1 .34 fim (H2OA), 
1.79 fim (H2OC), 1.96 ftm (H2OD) and at 1.2 fim from the FeH 
band) were derived for TWA22 A and B (see Fig.|9]l. The results 
were compared to the values computed from the ML03 and G05 
spectral libraries of field dwarfs. They were also compared to 
values derived for young dwarfs (Slesnick et al.'(2004); Lodieu 
|et al.[ ( [2008| l, hereafter S04 and L08) to test the sensitivity of 
these indexes to surface gravity (age). In fact. The H2O D and 
FeH indexes values tend to increase with age for M5-L2 dwarfs, 
and could disturb our analysis. We then used a mean weight of 
the individual spectral type estimations from H2OA, H2OC and 



the recently defined AUers H2O index at 1.55 fim (see AUers 
[et al.[ ( [2007| )) to infer M5 ± 1 and M5.1 ± 1 specti-al types for 
TWA22 A and B respectively. These results mach the M6 ±1.5 
and M5 + 1 values derived for TWA22 A and B from the H2OD 
and FeH indexes for the 2 objects. Based on the /T-band photom- 
etry and the associated bolometric corrections of j^Golimowski 



et al 



( [2004) l, we derive a luminosity of tog('L/Loj=-2.11+0.13 
dexTor TWA22 A and /o§('//Lo)=-2.30±0.16 dex TWA22 B 
Using the T(,yy-spectral type conversion scales for intermediate 



gravity objects ( Luhman et al.|2003 l, we find an initial estima- 
tion of Teff - 2990;^jgQ K for both components. 
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Fig. 5. Same as Fig. [ijbut for H band. In this case, differences be- 
tween M6V field dwarf spectrum and TWA22 spectra are impor- 
tant. This could arise from low gravity or flux losses introduced 
either in the standard star datacube and during the extraction pro- 
cess. 



4.3.3. Study of narrow lines 

Depths of many narrow lines were studied to provide additional 
information on the surface gravity of TWA22AB and particu- 
larly on its age. Following the Sembach & Savage ( 1992| l method 
to measure pseudo-equivalent widths and their associated uncer- 
tainties, we derived the equivalent widths of strong atomic lines 
over the J, H and K bands. They were computed for narrow 
lines at 1.106 /zm (TiO and H2O), 1.220 yum (FeH - Q branch), 
1.313 /im (Al I) and 1.626 yum (OH), and for the K I doublets at 
1.169, 1.177, 1.243 and 1.253 /vm. The results were compared 
with pseudo-equivalent widths of old field dwarfs (COS, ML03) 
and young Upper Sco dwarfs (S04, LOS). The use of both li- 
brairies confirmed the strong surface gravity dependency of the 
K I lines, more moderate for the Al 1, FeH and OH lines. Due 
to the degeneracy in terms of effective temperature and surface 
gravity, pseudo-equivalent widths alone are not sufficient for a 
precise spectral type determination of TWA22 AB. They remain 
however compatible with narrow lines depths of young and old 
dwarfs of spectral types later than M4. 

If we now assume a spectral type M6 ± 1 for both compo- 
nents, the pseudo-equivalent widths of TWA22 AB appear inter- 
mediate between values found for young and field dwarfs (see 
Fig. [TO) l. This is confirmed using a visual comparison with an 
evolutionary sequence of M6 dwarfs composed of the old field 
dwarf GL 406, the intermediate old companion AB Doc C (Age 
75 Myr, M5.5, [Close et alT] ( |2007| l) and a young M5.5 dwarf 
from the Orion nebulae (age ~ 1-2 Myr, S04). Together with the 
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Fig. 6. Same as Fig. |4]but for the K band. Our spectra are still 
contaminated by 50 Hz pick-up noise residuals. This effects is 
noticible around 2.3 ^m. Spectra look like M6V field dwarf. 



other age indicators, these intermediate surface gravity features 
confirm that TWA22 AB is likely to be a young binary system. 
However, their uncertainties remain significantly large to not as- 
sign a precise age. 



4.4. Gravity and effective temperature from atmospheric 
models 

For a fine determination of the effective temperatures and sur- 
face gravities of TWA22 A and B, we compared our observed 
spectra with theoretical templates from the GALA model v2.6.1 
(Brott & Hauschildt 2005). This library is updated from I Allard 
et al. (2001 ). It benefits from improved molecular dissociation 
constants, additional dust species with opacities, spherical sym- 
metry, and a mixing length parameter 2.0 x Hp. The temperature 
ranges in the templates from 2000 to 10 000 K and the gravity 
from -0.5 to 5.5 but we limited our analysis to 2000^r < T^ff < 
4000A' and 3.5 < log(g) < 5.5. Theoretical spectra were con- 
volved with a gaussian to match the SlNFONl spectral resolu- 
tion and interpolated to the TWA22 AB wavelength grid. Least 
squares minimization was applied to find templates that fit the 
TWA22 A and B continuum avoiding zones polluted by remain- 
ing oscillations. 

The TWA22 A least-square map in the J band constrains 
the temperature between 2800 to 3100 K and is minimized for 
log(g)=4.5 and Tp/y =2900 K. In H+K band, our minimization 
failed to reproduce faithfuly the TWA22 A spectra and makes us 
suspect the existence of a constant flux loss in H band during 
the spectral extraction process. To limit this systematic effect. 
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Fig. 7. Comparison of the TWA22 B J band spectrum (red) to 
spectra (black) of young Upper Sco and Orion nebulae cluster 
objects ( |Slesnick et al.||2004[ [Lodieu et aL||2008| ). We clearly 
notice that the TWA22 B spectral slope is redder than that of ref- 
erence spectra. The TWA22 B spectrum is very similar to those 
of young M6 and M7 dwarfs. We reported identified atomic fea- 
tures (blue). Molecular absorptions (FeH bands were identified 
by Gushing et al.| (i2QQ3 )) are flagged in green and telluric resid- 
uals in orange. Atomic absorptions are indicated in blue. 




2.2 

Wavelengh (yum) 

Fig. 8. Comparison of the TWA22 B K band spectrum (red) 
to spectra of young Upper Sco standards (black) at R=350 
( Gorlova et al. 2003 1. The TWA22 B spectrum was convolved 
with a gaussian to match the resolution of standard star spectra. 
Our spectrum is reproduced by spectra of young M5.5 and M7 
dwarf. 




Fig. 9. H2O spectral indexes computed on libraries of field 
( McLean et al.| 2003[ |Cushing et al. 2005)1 a nd young dwarfs 
spectra (Slesnick et al. 2004| [Lodieu et al.[ 2008 ). Redundancies 
between libraries have been checked. Young dwarfs M4-L2 
dwarfs follow trends of field dwarfs exept for H2O D. The dis- 
persion of values for young upper Sco brown dwarfs does not 
seem to arise from redenning effects. When overplotted, TWA22 
A and B values point an M6 ± 1 spectral type. In the upper right 
comer of each plot are reported the Icr errors derived from the 
residuals to the linear fit. 



the minimization was applied separately in H and K bands. In 
//-band, the effective temperature is minimized between 2600 K 
to 3000 K in the full space of surface gravities explored. The 
minimum is located at 2800 K and log(g)=4.5. The K band is 
well reproduced by 2900 and 3000 K templates irrespective of 
gravity. Summing the three bands, we estimate an effective tem- 
perature Tf// =2900+^11° K for TWA22 A. Conducting a similar 
analysis for the component TWA22 B, we derive an effective 



Luhman et al. 



(2003 I 



temperature r,.//=2900 K. Using the^ 

scale, these temperatures estimations correspond respectively to 
M7+^ and M7+°^ spectral types for TWA22 A and B. This is also 
consistent with spectral types estimated in part. 4.3.2. 

For a fine determination of the surface gravity from synthetic 
spectra, we computed the equivalent widths of K I lines in the J 
band on each spectra of TWA22 A and B. We then compared 
the values to TWA22 A and B to restrain the acceptable gravity 
domain (see Fig. [TT] i. We then estimate that the surface gravity 
is located between log(g)=4.0 and 5.5 for TWA22 A and B. 
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Fig. 10. Equivalent widths in the K I lines computed on libraries 
of young and field dwarfs spectra. Young dwarfs spectra have 
weak K I lines and therefore low equivalent widths compared 
to field dwarfs as a consequence of their low surface gravity. 
TWA22 A and B values are reported as red and green lines re- 
spectively with their associated uncertainties (dashed lines). K 
I equivalent widths of field dwarfs point a M5V spectral type 
different from our M6 ± 1 spectral type estimation from contin- 
uum fitting. This effect is probably the result of weakening of K 
I lines associated to the estimated very young age of the system. 



5. Discussion 

5. 1 . Evolutionary models predictions 

The membership of TWA22 AB to TW Hydrae constrains the 



age of the system to 3-20 Myr (Barrado Y Navascues 2006 



[Scholz et al.|2007t |de la Reza et al.|2006| l. Based on our astro- 
metric observations combined with an accurate distance deter- 
mination, we were able to derive the dynamical mass of this tight 
binary. From photometry and spectroscopy, we derived near- 
IR fluxes, luminosity, spectral type, effective temperatures and 
the surface gravity of each component. Finally, spectroscopy 
tends to indicate that both components have intermediate sur- 
face gravity features in their spectra, supporting a young age for 
TWA22 AB. Assuming the TWA age for this system, we can now 
compare the measured total dynamical mass of the binary with 




Log(g) 



Fig. 11. Iso-contours plots of K I lines equivalent widths com- 
puted on each spectral templates of the GAIA library v2.6.1 Red 
color indicates high values. The contours for TWA22 A (red) 
and B (blue) pseudo-equivalents widths values are overplotted. 
The long-dashed lines represent limits on TWA22 A (red) and 
TWA22 B (blue) temperatures. Gravity is estimated inside these 
temperatures boxes. 



the total mass predicted by evolutionary models of Baraffe et al. 
(1998; hereafter BCAH98). Model predictions are based on the 
JHK photometry, the luminosity and the effective temperature 
of both components (see Fig [H]!. At 8 Myr, BCAH98 models 
systematically under-estimate the total mass by a factor of ~ 2. 
This factor varies from 3 to 1.3 between 3-5 Myr and 20 Myr. 
The mass is still strongly under-estimated using other evolution- 



ary models of very low mass stars (D'Antona & Mazzitelli 1994 



1997[ |. Alternatively, if we artificially change the system age to 
30 Myr, the model predictions match relatively well our obser- 
vations. 

The apparent discrepancy between our observations and the 
model predictions at the age of TWA leads us to consider four 
explanations: 

1. Remaining uncertainties are present in our data reduction 
and interpretation related to the astrometry, photometry and 
spectroscopy extraction process, the atmosphere model used 
or the assumption on the system itself, 

2. The system is of higher multiplicity than observed. 
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Fig. 12. Top-Left : Confrontation of the binary direct mass mea- 
surement to predictions of the BCAH98 masses for different 
ages from Mk- Different estimations of the age of TW Hydrae 
(|de la Reza et al.|2006HBmado Y Navascues 20061 [Scholz et 
2007 1 are reported on the graph. Errors on the photometry are 
propagated on predictions (dotted lines). Bottom-Left : Same 
as top-left but for predictions from Mh- Top-right : Same as 
Top-Left but for predictions from Mj. Bottom-right : Same as 
Top-Left but for predictions from our estimated TWA22 A/B lu- 
minosities. In this plan, predictions from BCAH98, DM94 and 
DM97 models are nearly the same within our uncertainties. 

3. Evolutionary model predictions are correct and the age es- 
timate of TWA22AB is cuiTendy incoiTect. TWA22AB 
would be then slightly older and aged of 30 Myr, 

4. Finally, the TWA22 AB age is 8 Myr and evolutionary mod- 
els themselves do not predict correctly the physical proper- 
ties of very low mass stars at young ages. 

Before drawing important conclusions on the validity of evo- 
lutionary models at young ages and very low masses, we con- 
sider below the three first explanations. 

5.2. Data reduction and interpretation uncertainties 

Systematics on the estimation of the relative position and near- 
IR fluxes of TWA22 A and B seems very unlikely. Our analysis 
relies on the use of several imaging analysis techniques (aper- 
ture photometry, PSF fitting, deconvolution), already used and 
tested in various contexts. The tight binary TWA22 AB does not 
represent itself a difficult case. In addition, at each epoch, consis- 
tent results were found on several observing sequences obtained 
during the night. 

Systematics in the spectroscopic observation and extraction 
seem more probable for the determination of effective temper- 
ature and surface gravity. Differential flux losses over the J or 
H + K spectral range may have occurred due to the limited size 
of the SINFONI field of view. The impact of this effect can be 
simulated by adding a linear slope in our spectral minimization 
over the different spectral bands. The results do not change sig- 



nificantly our analysis based on continuum fitting or spectral 
indexes. It does not affect at all the study of narrow lines and 
our surface gravity estimation. A non-linear differential flux loss 
could be responsible for our failure to faithfully reproduce the 
TWA22 A and B spectra in H-band using either empirical or syn- 
thetic libraries. Finally, the atmosphere models were also used in 
various conditions (metallicity, mixing length, different opacity 
tables) without drastically changing our results. 



5.3. Higher multiplicity hypothesis 

Considering our data reduction and analysis as robust, we may 
wonder whether our basic assumptions concerning the system 
itself are correct. Actually, we cannot exclude from our observa- 
tions that TWA22 AB is a multiple system of higher order. One 
or even both components could be in fact unresolved binaries. 
In such case, the derived effective temperatures as well as the 
estimated spectral types would not be strongly modified. 

Dynamically speaking, the stability of the system would re- 
quire the separation of the individual sub-components to be sig- 
nificantly less than the size of the main orbit, typically by a factor 
3-4 dA rtymowicz & Lubow|199 4). To refine this estimate in the 
present case, we performed 3-body simulations using the sym- 
plectic code HJS (| Beust|2003 l) dedicated to hierarchical systems. 
We assume the fitted orbit and split one of the two components 
into 2 equal mass bodies, with a coplanar orbit with repect to 
the wide orbit and a given semi-major axis, and assuming initial 
zero eccentricity. We find that the system remains stable up to a 
separation of ~ 0.4 AU between the two sub-components. Note 
that this is a priori the most stable configuration, as a split into 
unequal masses would lead to less stability for the lighter com- 
ponent. We also checked that highly inclined relative configu- 
rations are physically unstable: they inevitably lead to a strong 
Kozai resonance characteristic for such triple systems ( |Kozai| 



1962 Beust et al. 1997 1 that cause the eccentricity of the in- 
ner orbit to be pumped up to ~ 1, leading to a physical colli- 
sion between the two individual components. Hence 0.4 AU can 
be considered as the widest possible separation for hypotheti- 



cal sub-components. This is in agreement with Artymowicz & 
Lubow (1994). 

0.4 AU (-22 mas) remains below the PSF of the VLT/NACO 
images (~ 1 AU given the distance of TWA22). This would ex- 
plain why the inner binary would not be resolved. However, we 
did not notice any PSF-lengthening in the images. With a sepa- 
ration of 0.4 AU we could expect to see one. Does it suggest that 
the actual separation is significantly smaller ? 

An orbit with 0.4 AU separation would correspond to an or- 
bital period of 0.8 yr and a radial velocity wobble of ±7 km s"' if 
we take into account the 27° inclination with respect to the plane 
of the sky. Even unlikely, this modulation could not to have been 
detected during the monitoring (split up into two periods of 6 
and 1 months). But if we assume a separation of ~ 0.1 AU to be 
compatible with the absence of PSF-lengthening, now the radial 
velocity wobble jumps to +14kms"' over a 0. 1 yr period. Such 
a variation was not detected in the radial velocity dataset (see 
I Teixeira et al.| ( 2009] l). Finally, no photocenter scatter is present 
around our two-body orbital solution. A motion of ~5 mas is 
expected along the orbital period for a separation of ~ 0. 1 AU 
between the subcomponents. 

Ultimately, we cannot definitely rule out the possibility that 
at least one of the two components of TWA22 is itself a binary, 
consisting of two nearly equal mass bodies. But combined dy- 
namical and observational constraints show that the range of 
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possible separations is fairly narrow, typically 0.1-0.2 AU. Also 
the system needs to be at least roughly coplanar. 

5.4. Age and membership of TWA22AB 

Given the good agreement between observations and model pre- 
dictions at 30 Myr, we can consider that the cuiTent age esti- 
mate of TWA22 AB is possibly incorrect. This age is currently 
infered from the membership to TWA. Since the age of TWA is 
well established at 8 Myr from various age diagnostics, a reliable 
explanation concerns the membership to TWA itself. 

S03 identified TWA22 as a new member of TWA mainly 
from the observed Li absorption line at 6708 A and Ha emis- 
sion line. The Li line is stronger (EW=510 mA) than those of 
early-M dwarfs members of Pic and leads S03 to suggest an 
age <10 Myr (see Fig. 8 of S03). They derived in addition a 
photometric distance of 22 pc for TWA22, confirming the prox- 
imity of this young system. More recently, Mamajek (2005) dis- 
cussed the membership of TWA22 AB to TWA based on its kine- 
matics properties. Mamajek (2005) estimates a probability of 
2% for TWA22 to be a member of TWA from an implemented 
convergence point technique ( de Bruijne|1999| . However, Song 
[eTaLl ( p006| mentioned that the strong Li line of TWA22 AB 
is observed only for young active M dwarfs in the direction of 
TWA with the exception of a very few M-type members of the /3 
Pictoris moving group (BPMG). Finally, Mentuch et al.| ( |200§] l 
obtained a new visible high resolution spectrum of TWA22 AB. 
They confirmed the strong equivalent line of the 6708 A lithium 
absorption (EW=616 ±21 mA; the strongest measured in their 
sample composed of young association members). They esti- 
mated Teff of 2990 + 13 K (compatible with the individual Teff 
derived in Part 4.4) and log(g)=4.2 ± 0.05 dex for the unresolved 
system. These new elements tend to confirm that TWA22 is a 
young system (age < 30 Myr; see BCAH98 predictions). 

To conciliate past and present results, we can consider the 
possibility that TWA22 AB is a member of the BPMG. With an 
M6 + 1 spectral type, TWA22 AB is probably close to the Li- 
depletion boundary (LDB) of TWA or/? Pic, which could possi- 
bly explain a significantly stronger EW(Li) than those observed 
for early-M dwarfs of these two associations. We can also no- 
tice that the A 6708 A line shows some variations between the 
S03 and M08 measurments. In addition, the observed EW(Hq,) of 
TWA22, used as a second indicator of youth, is compatible with 
those of GJ799 A and B, M4.5 members of /3 Pic Payawardhana| 
let al.|20061 l. 

Finally, the projected position of TWA22 AB reveals that the 
system is isolated from other members of TWA. Its distance is 
more compatible with the mean distance of the BPMG mem- 



bers. Teixeira et al. ( 2009 1 have recently measured the proper 
motion, the trigonometric parallax and the mean radial veloc- 
ity of TWA22 AB. They determined for the first time the helio- 
centric space motion of TWA22 AB. From a detailed kinematic 
analysis they did not ruled out TWA22 from TW Hydrae but they 
demonstrated that it was a more probable member of the BPMG. 



6. Conclusions 

NACO resolved for the first time the young object TWA 22 as 
a tight binary with a projected separation of 1.76 AU. 80 % of 
the binary orbit was covered duiing a 4 years observation pro- 
gram conducted with this instrument. We inferred a 220 ±21 
Mjiip total mass for the system and we obtained the individual 
magnitudes of each component in the near infrared. This places 



TWA22 A and B at the substellar boundary. We complete the 
characterization of the system components with medium res- 
olution individual SINFONI spectra in the J, H and K bands. 
Our spectra were compared with empirical library of young and 
field M dwarfs. We derived a M6 ± 1 spectral type from con- 
tinuum fitting, spectral indexes and equivalent widths. Spectral 
templates were also used to estimate Teff-29QQ + 200 K for 
TWA22 A and r,//=2900 k for TWA22 B while the suiface 
gravity was constrained to 4.0 < log(g) < 5.5 dex. These funda- 
mental properties can be directly compared with commonly used 
evolutionary tracks provided that the age of the system is known 
accurately. 

The age of TWA22 was still a mater of debate at the begin- 
ning of our study. TWA22 was reported as member of the young 
association TW Hydrae, and alternatively as a possible member 
of the BPMG. At the age of TW Hydrae and BPMG, the new 
and precious benchmark brought by this system seems to point 
an underestimation of the predicted mass from our photometry. 
However, the dynamical mass appears correctly estimated by the 
models if we consider a 30 Myr old system. This led us recon- 
sider the membership of TWA22. 

While the spectroscopy tends to confirm the youth of this 
system, a recent kinematic study rejected TWA22 as a member 
of the TW Hydrae and of the 30 Myr old Tucana-Horologium 
associations. It didn't exclude the membership of TWA22 to the 
BPMG. Also, we can not rule out the possibility that TWA22 is 
not associated with any of these moving groups. 

Finally, we does not exclude firmly that TWA22AB com- 
ponent could be in fact unresolved binaries with coplanar inner 
orbits characterized by semi-major axis lower than 0.4 AU. The 
models predictions would match the measured dynamical mass 
of a triple or quadruple system. In this context, future monitor- 
ing of TWA22AB with improved angular resolution could allow 
the resolution of the hypothetical inner binaries. 
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